3.5.98 \(\int \frac {x^3 (a+b \text {arccosh}(c x))}{(d+e x^2)^2} \, dx\) [498]

3.5.98.1 Optimal result
3.5.98.2 Mathematica [C] (verified)
3.5.98.3 Rubi [A] (verified)
3.5.98.4 Maple [C] (warning: unable to verify)
3.5.98.5 Fricas [F]
3.5.98.6 Sympy [F]
3.5.98.7 Maxima [F]
3.5.98.8 Giac [F(-2)]
3.5.98.9 Mupad [F(-1)]

3.5.98.1 Optimal result

Integrand size = 21, antiderivative size = 562 \[ \int \frac {x^3 (a+b \text {arccosh}(c x))}{\left (d+e x^2\right )^2} \, dx=\frac {d (a+b \text {arccosh}(c x))}{2 e^2 \left (d+e x^2\right )}-\frac {(a+b \text {arccosh}(c x))^2}{2 b e^2}-\frac {b c \sqrt {d} \sqrt {-1+c^2 x^2} \text {arctanh}\left (\frac {\sqrt {c^2 d+e} x}{\sqrt {d} \sqrt {-1+c^2 x^2}}\right )}{2 e^2 \sqrt {c^2 d+e} \sqrt {-1+c x} \sqrt {1+c x}}+\frac {(a+b \text {arccosh}(c x)) \log \left (1-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}-\sqrt {-c^2 d-e}}\right )}{2 e^2}+\frac {(a+b \text {arccosh}(c x)) \log \left (1+\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}-\sqrt {-c^2 d-e}}\right )}{2 e^2}+\frac {(a+b \text {arccosh}(c x)) \log \left (1-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}+\sqrt {-c^2 d-e}}\right )}{2 e^2}+\frac {(a+b \text {arccosh}(c x)) \log \left (1+\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}+\sqrt {-c^2 d-e}}\right )}{2 e^2}+\frac {b \operatorname {PolyLog}\left (2,-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}-\sqrt {-c^2 d-e}}\right )}{2 e^2}+\frac {b \operatorname {PolyLog}\left (2,\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}-\sqrt {-c^2 d-e}}\right )}{2 e^2}+\frac {b \operatorname {PolyLog}\left (2,-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}+\sqrt {-c^2 d-e}}\right )}{2 e^2}+\frac {b \operatorname {PolyLog}\left (2,\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}+\sqrt {-c^2 d-e}}\right )}{2 e^2} \]

output
1/2*d*(a+b*arccosh(c*x))/e^2/(e*x^2+d)-1/2*(a+b*arccosh(c*x))^2/b/e^2+1/2* 
(a+b*arccosh(c*x))*ln(1-(c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))*e^(1/2)/(c*(-d)^ 
(1/2)-(-c^2*d-e)^(1/2)))/e^2+1/2*(a+b*arccosh(c*x))*ln(1+(c*x+(c*x-1)^(1/2 
)*(c*x+1)^(1/2))*e^(1/2)/(c*(-d)^(1/2)-(-c^2*d-e)^(1/2)))/e^2+1/2*(a+b*arc 
cosh(c*x))*ln(1-(c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))*e^(1/2)/(c*(-d)^(1/2)+(- 
c^2*d-e)^(1/2)))/e^2+1/2*(a+b*arccosh(c*x))*ln(1+(c*x+(c*x-1)^(1/2)*(c*x+1 
)^(1/2))*e^(1/2)/(c*(-d)^(1/2)+(-c^2*d-e)^(1/2)))/e^2+1/2*b*polylog(2,-(c* 
x+(c*x-1)^(1/2)*(c*x+1)^(1/2))*e^(1/2)/(c*(-d)^(1/2)-(-c^2*d-e)^(1/2)))/e^ 
2+1/2*b*polylog(2,(c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))*e^(1/2)/(c*(-d)^(1/2)- 
(-c^2*d-e)^(1/2)))/e^2+1/2*b*polylog(2,-(c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))* 
e^(1/2)/(c*(-d)^(1/2)+(-c^2*d-e)^(1/2)))/e^2+1/2*b*polylog(2,(c*x+(c*x-1)^ 
(1/2)*(c*x+1)^(1/2))*e^(1/2)/(c*(-d)^(1/2)+(-c^2*d-e)^(1/2)))/e^2-1/2*b*c* 
arctanh(x*(c^2*d+e)^(1/2)/d^(1/2)/(c^2*x^2-1)^(1/2))*d^(1/2)*(c^2*x^2-1)^( 
1/2)/e^2/(c^2*d+e)^(1/2)/(c*x-1)^(1/2)/(c*x+1)^(1/2)
 
3.5.98.2 Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 1.38 (sec) , antiderivative size = 694, normalized size of antiderivative = 1.23 \[ \int \frac {x^3 (a+b \text {arccosh}(c x))}{\left (d+e x^2\right )^2} \, dx=\frac {\frac {2 a d}{d+e x^2}+2 a \log \left (d+e x^2\right )+b \left (-2 \text {arccosh}(c x)^2+2 \text {arccosh}(c x) \left (\log \left (1+\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{-i c \sqrt {d}+\sqrt {-c^2 d-e}}\right )+\log \left (1-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{i c \sqrt {d}+\sqrt {-c^2 d-e}}\right )\right )+2 \text {arccosh}(c x) \left (\log \left (1+\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{i c \sqrt {d}-\sqrt {-c^2 d-e}}\right )+\log \left (1+\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{i c \sqrt {d}+\sqrt {-c^2 d-e}}\right )\right )-i \sqrt {d} \left (\frac {\text {arccosh}(c x)}{-i \sqrt {d}+\sqrt {e} x}+\frac {c \log \left (\frac {2 e \left (i \sqrt {e}+c^2 \sqrt {d} x-i \sqrt {-c^2 d-e} \sqrt {-1+c x} \sqrt {1+c x}\right )}{c \sqrt {-c^2 d-e} \left (\sqrt {d}+i \sqrt {e} x\right )}\right )}{\sqrt {-c^2 d-e}}\right )-i \sqrt {d} \left (-\frac {\text {arccosh}(c x)}{i \sqrt {d}+\sqrt {e} x}-\frac {c \log \left (\frac {2 e \left (-\sqrt {e}-i c^2 \sqrt {d} x+\sqrt {-c^2 d-e} \sqrt {-1+c x} \sqrt {1+c x}\right )}{c \sqrt {-c^2 d-e} \left (i \sqrt {d}+\sqrt {e} x\right )}\right )}{\sqrt {-c^2 d-e}}\right )+2 \operatorname {PolyLog}\left (2,\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{i c \sqrt {d}-\sqrt {-c^2 d-e}}\right )+2 \operatorname {PolyLog}\left (2,\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{-i c \sqrt {d}+\sqrt {-c^2 d-e}}\right )+2 \operatorname {PolyLog}\left (2,-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{i c \sqrt {d}+\sqrt {-c^2 d-e}}\right )+2 \operatorname {PolyLog}\left (2,\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{i c \sqrt {d}+\sqrt {-c^2 d-e}}\right )\right )}{4 e^2} \]

input
Integrate[(x^3*(a + b*ArcCosh[c*x]))/(d + e*x^2)^2,x]
 
output
((2*a*d)/(d + e*x^2) + 2*a*Log[d + e*x^2] + b*(-2*ArcCosh[c*x]^2 + 2*ArcCo 
sh[c*x]*(Log[1 + (Sqrt[e]*E^ArcCosh[c*x])/((-I)*c*Sqrt[d] + Sqrt[-(c^2*d) 
- e])] + Log[1 - (Sqrt[e]*E^ArcCosh[c*x])/(I*c*Sqrt[d] + Sqrt[-(c^2*d) - e 
])]) + 2*ArcCosh[c*x]*(Log[1 + (Sqrt[e]*E^ArcCosh[c*x])/(I*c*Sqrt[d] - Sqr 
t[-(c^2*d) - e])] + Log[1 + (Sqrt[e]*E^ArcCosh[c*x])/(I*c*Sqrt[d] + Sqrt[- 
(c^2*d) - e])]) - I*Sqrt[d]*(ArcCosh[c*x]/((-I)*Sqrt[d] + Sqrt[e]*x) + (c* 
Log[(2*e*(I*Sqrt[e] + c^2*Sqrt[d]*x - I*Sqrt[-(c^2*d) - e]*Sqrt[-1 + c*x]* 
Sqrt[1 + c*x]))/(c*Sqrt[-(c^2*d) - e]*(Sqrt[d] + I*Sqrt[e]*x))])/Sqrt[-(c^ 
2*d) - e]) - I*Sqrt[d]*(-(ArcCosh[c*x]/(I*Sqrt[d] + Sqrt[e]*x)) - (c*Log[( 
2*e*(-Sqrt[e] - I*c^2*Sqrt[d]*x + Sqrt[-(c^2*d) - e]*Sqrt[-1 + c*x]*Sqrt[1 
 + c*x]))/(c*Sqrt[-(c^2*d) - e]*(I*Sqrt[d] + Sqrt[e]*x))])/Sqrt[-(c^2*d) - 
 e]) + 2*PolyLog[2, (Sqrt[e]*E^ArcCosh[c*x])/(I*c*Sqrt[d] - Sqrt[-(c^2*d) 
- e])] + 2*PolyLog[2, (Sqrt[e]*E^ArcCosh[c*x])/((-I)*c*Sqrt[d] + Sqrt[-(c^ 
2*d) - e])] + 2*PolyLog[2, -((Sqrt[e]*E^ArcCosh[c*x])/(I*c*Sqrt[d] + Sqrt[ 
-(c^2*d) - e]))] + 2*PolyLog[2, (Sqrt[e]*E^ArcCosh[c*x])/(I*c*Sqrt[d] + Sq 
rt[-(c^2*d) - e])]))/(4*e^2)
 
3.5.98.3 Rubi [A] (verified)

Time = 1.39 (sec) , antiderivative size = 562, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.095, Rules used = {6374, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {x^3 (a+b \text {arccosh}(c x))}{\left (d+e x^2\right )^2} \, dx\)

\(\Big \downarrow \) 6374

\(\displaystyle \int \left (\frac {x (a+b \text {arccosh}(c x))}{e \left (d+e x^2\right )}-\frac {d x (a+b \text {arccosh}(c x))}{e \left (d+e x^2\right )^2}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {(a+b \text {arccosh}(c x)) \log \left (1-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}-\sqrt {c^2 (-d)-e}}\right )}{2 e^2}+\frac {(a+b \text {arccosh}(c x)) \log \left (\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}-\sqrt {c^2 (-d)-e}}+1\right )}{2 e^2}+\frac {(a+b \text {arccosh}(c x)) \log \left (1-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{\sqrt {c^2 (-d)-e}+c \sqrt {-d}}\right )}{2 e^2}+\frac {(a+b \text {arccosh}(c x)) \log \left (\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{\sqrt {c^2 (-d)-e}+c \sqrt {-d}}+1\right )}{2 e^2}+\frac {d (a+b \text {arccosh}(c x))}{2 e^2 \left (d+e x^2\right )}-\frac {(a+b \text {arccosh}(c x))^2}{2 b e^2}+\frac {b \operatorname {PolyLog}\left (2,-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}-\sqrt {-d c^2-e}}\right )}{2 e^2}+\frac {b \operatorname {PolyLog}\left (2,\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{c \sqrt {-d}-\sqrt {-d c^2-e}}\right )}{2 e^2}+\frac {b \operatorname {PolyLog}\left (2,-\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{\sqrt {-d} c+\sqrt {-d c^2-e}}\right )}{2 e^2}+\frac {b \operatorname {PolyLog}\left (2,\frac {\sqrt {e} e^{\text {arccosh}(c x)}}{\sqrt {-d} c+\sqrt {-d c^2-e}}\right )}{2 e^2}-\frac {b c \sqrt {d} \sqrt {c^2 x^2-1} \text {arctanh}\left (\frac {x \sqrt {c^2 d+e}}{\sqrt {d} \sqrt {c^2 x^2-1}}\right )}{2 e^2 \sqrt {c x-1} \sqrt {c x+1} \sqrt {c^2 d+e}}\)

input
Int[(x^3*(a + b*ArcCosh[c*x]))/(d + e*x^2)^2,x]
 
output
(d*(a + b*ArcCosh[c*x]))/(2*e^2*(d + e*x^2)) - (a + b*ArcCosh[c*x])^2/(2*b 
*e^2) - (b*c*Sqrt[d]*Sqrt[-1 + c^2*x^2]*ArcTanh[(Sqrt[c^2*d + e]*x)/(Sqrt[ 
d]*Sqrt[-1 + c^2*x^2])])/(2*e^2*Sqrt[c^2*d + e]*Sqrt[-1 + c*x]*Sqrt[1 + c* 
x]) + ((a + b*ArcCosh[c*x])*Log[1 - (Sqrt[e]*E^ArcCosh[c*x])/(c*Sqrt[-d] - 
 Sqrt[-(c^2*d) - e])])/(2*e^2) + ((a + b*ArcCosh[c*x])*Log[1 + (Sqrt[e]*E^ 
ArcCosh[c*x])/(c*Sqrt[-d] - Sqrt[-(c^2*d) - e])])/(2*e^2) + ((a + b*ArcCos 
h[c*x])*Log[1 - (Sqrt[e]*E^ArcCosh[c*x])/(c*Sqrt[-d] + Sqrt[-(c^2*d) - e]) 
])/(2*e^2) + ((a + b*ArcCosh[c*x])*Log[1 + (Sqrt[e]*E^ArcCosh[c*x])/(c*Sqr 
t[-d] + Sqrt[-(c^2*d) - e])])/(2*e^2) + (b*PolyLog[2, -((Sqrt[e]*E^ArcCosh 
[c*x])/(c*Sqrt[-d] - Sqrt[-(c^2*d) - e]))])/(2*e^2) + (b*PolyLog[2, (Sqrt[ 
e]*E^ArcCosh[c*x])/(c*Sqrt[-d] - Sqrt[-(c^2*d) - e])])/(2*e^2) + (b*PolyLo 
g[2, -((Sqrt[e]*E^ArcCosh[c*x])/(c*Sqrt[-d] + Sqrt[-(c^2*d) - e]))])/(2*e^ 
2) + (b*PolyLog[2, (Sqrt[e]*E^ArcCosh[c*x])/(c*Sqrt[-d] + Sqrt[-(c^2*d) - 
e])])/(2*e^2)
 

3.5.98.3.1 Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 6374
Int[((a_.) + ArcCosh[(c_.)*(x_)]*(b_.))^(n_.)*((f_.)*(x_))^(m_.)*((d_) + (e 
_.)*(x_)^2)^(p_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*ArcCosh[c*x])^n, 
 (f*x)^m*(d + e*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[c^2*d 
 + e, 0] && IGtQ[n, 0] && IntegerQ[p] && IntegerQ[m]
 
3.5.98.4 Maple [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 3.01 (sec) , antiderivative size = 2132, normalized size of antiderivative = 3.79

method result size
derivativedivides \(\text {Expression too large to display}\) \(2132\)
default \(\text {Expression too large to display}\) \(2132\)
parts \(\text {Expression too large to display}\) \(2144\)

input
int(x^3*(a+b*arccosh(c*x))/(e*x^2+d)^2,x,method=_RETURNVERBOSE)
 
output
1/c^4*(1/2*a*c^6/e^2*d/(c^2*e*x^2+c^2*d)+1/2*a*c^4/e^2*ln(c^2*e*x^2+c^2*d) 
+b*c^4*(-1/8*(-2*(d*c^2*(c^2*d+e))^(1/2)*c^2*d+2*c^4*d^2+2*c^2*d*e-(d*c^2* 
(c^2*d+e))^(1/2)*e)/c^2/d/e^2/(c^2*d+e)*polylog(2,e*(c*x+(c*x-1)^(1/2)*(c* 
x+1)^(1/2))^2/(-2*c^2*d-2*(d*c^2*(c^2*d+e))^(1/2)-e))+(-2*(d*c^2*(c^2*d+e) 
)^(1/2)*c^2*d+2*c^4*d^2+2*c^2*d*e-(d*c^2*(c^2*d+e))^(1/2)*e)*c^2*d/e^4/(c^ 
2*d+e)*arccosh(c*x)^2-1/2*(-2*(d*c^2*(c^2*d+e))^(1/2)*c^2*d+2*c^4*d^2+2*c^ 
2*d*e-(d*c^2*(c^2*d+e))^(1/2)*e)*c^2*d/e^4/(c^2*d+e)*polylog(2,e*(c*x+(c*x 
-1)^(1/2)*(c*x+1)^(1/2))^2/(-2*c^2*d-2*(d*c^2*(c^2*d+e))^(1/2)-e))+1/4*(d* 
c^2*(c^2*d+e))^(1/2)/e/d/c^2/(c^2*d+e)*arccosh(c*x)^2+1/4*(-2*(d*c^2*(c^2* 
d+e))^(1/2)*c^2*d+2*c^4*d^2+2*c^2*d*e-(d*c^2*(c^2*d+e))^(1/2)*e)/c^2/d/e^2 
/(c^2*d+e)*arccosh(c*x)^2+(2*c^2*d-2*(d*c^2*(c^2*d+e))^(1/2)+e)/e^4*ln(1-e 
*(c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))^2/(-2*c^2*d-2*(d*c^2*(c^2*d+e))^(1/2)-e 
))*c^2*d*arccosh(c*x)+1/2*(d*c^2*(c^2*d+e))^(1/2)/e^2/(c^2*d+e)*arctanh(1/ 
4*(4*c^2*d+2*e*(c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))^2+2*e)/(c^4*d^2+c^2*d*e)^ 
(1/2))-1/4*(d*c^2*(c^2*d+e))^(1/2)/e^2/(c^2*d+e)*polylog(2,e*(c*x+(c*x-1)^ 
(1/2)*(c*x+1)^(1/2))^2/(-2*c^2*d+2*(d*c^2*(c^2*d+e))^(1/2)-e))-1/8*(d*c^2* 
(c^2*d+e))^(1/2)/e/d/c^2/(c^2*d+e)*polylog(2,e*(c*x+(c*x-1)^(1/2)*(c*x+1)^ 
(1/2))^2/(-2*c^2*d+2*(d*c^2*(c^2*d+e))^(1/2)-e))+(-2*(d*c^2*(c^2*d+e))^(1/ 
2)*c^2*d+2*c^4*d^2+2*c^2*d*e-(d*c^2*(c^2*d+e))^(1/2)*e)/e^3/(c^2*d+e)*arcc 
osh(c*x)^2+1/2*(2*c^2*d-2*(d*c^2*(c^2*d+e))^(1/2)+e)/e^3*ln(1-e*(c*x+(c...
 
3.5.98.5 Fricas [F]

\[ \int \frac {x^3 (a+b \text {arccosh}(c x))}{\left (d+e x^2\right )^2} \, dx=\int { \frac {{\left (b \operatorname {arcosh}\left (c x\right ) + a\right )} x^{3}}{{\left (e x^{2} + d\right )}^{2}} \,d x } \]

input
integrate(x^3*(a+b*arccosh(c*x))/(e*x^2+d)^2,x, algorithm="fricas")
 
output
integral((b*x^3*arccosh(c*x) + a*x^3)/(e^2*x^4 + 2*d*e*x^2 + d^2), x)
 
3.5.98.6 Sympy [F]

\[ \int \frac {x^3 (a+b \text {arccosh}(c x))}{\left (d+e x^2\right )^2} \, dx=\int \frac {x^{3} \left (a + b \operatorname {acosh}{\left (c x \right )}\right )}{\left (d + e x^{2}\right )^{2}}\, dx \]

input
integrate(x**3*(a+b*acosh(c*x))/(e*x**2+d)**2,x)
 
output
Integral(x**3*(a + b*acosh(c*x))/(d + e*x**2)**2, x)
 
3.5.98.7 Maxima [F]

\[ \int \frac {x^3 (a+b \text {arccosh}(c x))}{\left (d+e x^2\right )^2} \, dx=\int { \frac {{\left (b \operatorname {arcosh}\left (c x\right ) + a\right )} x^{3}}{{\left (e x^{2} + d\right )}^{2}} \,d x } \]

input
integrate(x^3*(a+b*arccosh(c*x))/(e*x^2+d)^2,x, algorithm="maxima")
 
output
1/2*a*(d/(e^3*x^2 + d*e^2) + log(e*x^2 + d)/e^2) + b*integrate(x^3*log(c*x 
 + sqrt(c*x + 1)*sqrt(c*x - 1))/(e^2*x^4 + 2*d*e*x^2 + d^2), x)
 
3.5.98.8 Giac [F(-2)]

Exception generated. \[ \int \frac {x^3 (a+b \text {arccosh}(c x))}{\left (d+e x^2\right )^2} \, dx=\text {Exception raised: TypeError} \]

input
integrate(x^3*(a+b*arccosh(c*x))/(e*x^2+d)^2,x, algorithm="giac")
 
output
Exception raised: TypeError >> an error occurred running a Giac command:IN 
PUT:sage2:=int(sage0,sageVARx):;OUTPUT:sym2poly/r2sym(const gen & e,const 
index_m & i,const vecteur & l) Error: Bad Argument Value
 
3.5.98.9 Mupad [F(-1)]

Timed out. \[ \int \frac {x^3 (a+b \text {arccosh}(c x))}{\left (d+e x^2\right )^2} \, dx=\int \frac {x^3\,\left (a+b\,\mathrm {acosh}\left (c\,x\right )\right )}{{\left (e\,x^2+d\right )}^2} \,d x \]

input
int((x^3*(a + b*acosh(c*x)))/(d + e*x^2)^2,x)
 
output
int((x^3*(a + b*acosh(c*x)))/(d + e*x^2)^2, x)